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I. INTRODUCTION 



In the last two decades two different thermodynamic theories have been applied in modeling fast non-equilibrium 
phenomena: the Extended Thermodynamics 0] and the Irreversible Thermodynamics with dynamic variables 
0,0. In the following investigations instead of referring directly to persons we sometimes distinct theories of Extended 
Thermodynamics (Rational Extended Thermodynamics , Extended Irreversible Thermodynamics 0, Q and Wave 
Approach of Thermodynamics 6,7]). 

Extended Thermodynamics (ET) lies on the kinetic theory of gases and postulates suitable balance equations for 
dissipative fluxes. In this way a hierarchical system in which the flux at step n becomes the wanted field at step n+1, 
may be obtained. Furthermore, the evolutionary systems are required to be symmetric hyperbolic so that finite speed 
of propagation of disturbances and well posedness of the Cauchy problem are guaranteed jlj . 

Irreversible Thermodynamics (IT) enlarges the basic state space by introducing certain additional dynamic variables 
whose evolution is ruled either by ordinary differential equations (local theory) or by partial differential equations 
(weakly nonlocal theory) 9, lOj . In such a framework the parabolic models are also allowed in that, once their solutions 
are interpreted in the light of the experimental measurements, they also lead to finite speed of propagation • ET 
is a special case of IT with internal variables, where the form of the evolution equations is restricted and the role of 
the internal variables is fixed. 

In the present paper we compare these two approaches by deriving conditions under which they lead to the same 
system of equations. To achieve that task we apply the classical Onsager analysis of Second Law P.[l^. 

We point out the central role of the entropy current in the two theories. In Rational Extended Thermodynamics 
the entropy current is a derived quantity. Its form is restricted by the balance structure of the evolution equations 
through the Second Law (see e.g. page 24 or page 60). The restriction can be expressed by the following differential 
expression: 

dj. = ^df. (I.l) 

Here ji is the ith component of the entropy current, s is the entropy density, a and j° denotes the densities their 
currents in the corresponding balances a -\- jf^ = 0, respectively. Here the entropy current appears as a potential 
on the space of the currents of the extensive quantities. Considering that the currents of the extensives j"' are 
also constitutive quantities, the above expression constraints considerably their form. The surprisingly restrictive 
constraints were proved to be compatible with some kinetic theory models. As a price some of the higher order 
currents were considered as constitutive quantities instead of being variables. 
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On the other hand in IT with dynamic variables |j] the starting point is the following form of the entropy current 

= (1.2) 

This form was suggested as a direct generalization of the classical heat flux over temperature expression ji = qi /T 
0- . 

The entropy current 1)1.1(1 is the consequence of the balance structure of the evolution equations. On the other hand 
from any given form of the entropy current one can determine, or at least restrict, the form of the evolution equations 
of the dynamic variables. In Extended Irreversible Thermodynamics one can meet both forms, J ou, Casas-Vasquez 
and Lebon seem to hesitate on the advantages and disadvantages of the two suggestions (compare [l^) Q P- 138 and 
[2, p. 56-58, [H). 

In the following we will investigate the compatibility of the Verhas form of the entropy current with the balance 
structure of the evolution equations in a general and highly nonlinear case. We give conditions when the balance 
structure and the HI.2|I entropy current are compatible. Our investigations show that if we do not insist strictly to 
the balance structure then the momentum hierarchy can be closed easily on the phenomenological level. 

In section 2, after considering a rigid body at rest and out of local equilibrium, we specify the basic balance 
equations together with the main constitutive assumptions. In section 3 wc obtain the evolution equations for the 
dynamic variables. In section 4 the conditions are pointed out under which these equations yield a 4 moments 
model and a 13 moments model of rigid heat conductor of Extended Thermodynamics. In section 5 we consider 
the propagation of heat waves at low temperatures. We reinspect some extended thermodynamic models 16] and 
determine the conditions under which they can be recovered in the framework of the presented theory. A concluding 
discussion of the obtained results is developed in section 6. 



II. BALANCE EQUATIONS AND SECOND LAW OF THERMODYNAMICS 

As a simple but representative example let us consider a rigid heat conductor at rest and let us start from the 
following local balance equation of the internal energy 

e + q^H^O, (II.l) 

where e is the density of internal energy, qi i = 1,2,3 are the components of the heat flux q, f = f>i = Xi 
i = 1, 2, 3 are the Cartesian coordinates of the points of the body and the Einstein convention of summation over the 
repeated indices has been applied. The local balance of entropy is given by 

s + ji'i^os, (II-2) 

with s standing for the entropy density, ji i = 1,2,3 for the components of the entropy current j and Us for the 
density of entropy production. Second Law of Thermodynamics forces as to be nonnegative. The only equilibrium 
variable will be the internal energy e, while the first dynamic variable is supposed to be the heat flux itself. Such 
an assumption, following by the celebrated Cattaneo's equation is included within the basic postulates of both 
theories under examination (see Chapter 1, and Q). Hence, it also will be our basic starting point. As a further 
dynamic variable let us choose a second order tensor # whose components will be denoted by , i,j — 1, 2, 3. The 
physical meaning of <I>y will be specified in section 5 on the example of second sound wave propagation. 
Hence we postulate the constitutive equation as 

F = F*(e,(7„$y), (H.3) 

where F means an element of the set of constitutive functions {s, j}. The basic thermodynamic state space is spanned 
by the 13 unknown quantities {e,(7i,$y } but later on some enlargements will be considered. Beside Eg. 1(11. 1(1 . the 
system of evolution equations for the 12 wanted fields qi and ^ij is needed. In order to derive such a system let us 
calculate as according to the classical procedures of Irreversible Thermodynamics. Along with Gyarmati and Verhas 
[^Il3j| we represent the entropy function, out of local equilibrium, as 

s(e, qi, $y ) So(e) + ^m.jqiqj + ^n.ijki^ij^ki- (H.4) 

In l|II.4|) so(e) is the local equilibrium entropy while the matrices of thermodynamic inductivities rriij and riijki are 
constitutive functions depending on the basic fields e, qi and ^ij . Moreover, we assume that the local equilibrium is 
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globally thermodynamically stable and there are no phase boundaries at the non-equilibrium part of the state space. 
Therefore, due to the principle of maximum entropy at the equilibrium, and riiji^i are negative definite. As the 
antisymmetric part of the inductivities m and n do not contribute to the entropy we assume that they are symmetric 
as follows: 

rriij = ruji, n^jki = n^uj- (11.5) 
Finally, along with Verhas [T^ IT3. ITsj l we represent the entropy current as 

ds ds ^ ^, 

Let us observe that the previous form of the entropy current is not the most general one and some extensions are 
also possible (see for instance ^3). Due to ljIL4|) . ji specialize to 

^5 1 1 

ji = "^^^ + -^Ahkqhqkqi + -^Shkim^hk^imqi + (n.7) 

^Chk]qhqk^ji+ + mhjqh^ji + ^Dhkhnj^hk^im^jz- 



and 



^ dnihk ^ dnhkim . ^ ^ dm^k „ ^ duhkir. 

oe oe oqj aqj 



In the next section we will use (|IL4ll and (|II.7|I in order to calculate the most important quantity related to this 
procedure, namely, the density of entropy production. Let us observe the special dependence of vector j on the 
components qi and of the basic field. Due to (jll.7|l the scalar jVi does not contain the divergence of any third 
order tensor to be interpreted as the flux of $. As far as the form of the evolution equations is concerned, such a 
property results to be crucial. 

III. BILINEAR FORM OF THE ENTROPY PRODUCTION 

Our aim here is to prove that the entropy production (Js can be put in the classical bilinear form of Irreversible 
Thermodynamics lU]. Some straightforward calculations show that, due to (jll.4p and (|II.7|I the function as can be 
written as follows 

1 . . 1 . 

— e'i + -A.kqi'iqk + Ahkqhqk'i + -j^Auk' iqhqk + 

--^Biiklm'i^hk^hn + -^{Bliklm + Bimhk)'^hk' i'^lm + -^Cikl'mqk^lTn + 
iiqh'm Im'jn 



^A^jq^e + ^Cijkqjqk + m[ik]qk + ^Fijki^kiq, 



-^Bhklm^lmqi'r 



+mhpqp'k + -^Dhklmp'i^lm^pi + Dhklmp^lTn'i^pi + 



+ -^^h.klmp^lm^pi'i + '^Bhklj'^lj^ + -^D h,kjlpqp^ jl + 
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hk- 



(IILl) 



In deriving (jlll.ip we have put 



^ ijkl 



drrii 



G 



dri; 



ijklpq 



ijkl 



pq 



A direct inspection of Ijlll.ip shows that each term stands for the combination of a relaxation and a transport. The 
quantities in the square brackets are regarded as process forces while their coefficients qi and $y as thermodynamic 
fluxes. Because the expression (jIII.I|l of the entropy production is of the usual bilinear form, let us apply Onsager's 
linear equations. In this way we get the evolution equations we are looking for, i. e. 



ihk 



LikQk + Miki^ki 

Phkiqi + Qhkij^i: 



(III.2) 
(III.3) 



where the symbols [q].; and [^]hk stand for the coefhcients of qi and ^hk respectively and the tensor functions L, 
M, P, Q are constitutive functions depending on the basic fields. With these conditions the above form is the solution 
of the entropy inequality. 

It is important to remark here, that the classical Coleman-Noll [23| and Liu |23| procedures for the exploitation of 
Second Law are not suitable to obtain the equations above since these techniques lead to thermodynamic restrictions 
on the constitutive equations which are not in the form of phenomenological evolution equations. After applying these 
exploitation techniques one can get a restricted form of the dissipation inequality. The Onsagerian approach could 
start here, with the identification of the fluxes and the forces. Indeed, the essence of the Onsager approach is the 
strict distinction between the undetermined constitutive quantities (currents) and given functions of the constitutive 
space (forces) . With a suitable distinction and some continuity requirements one can give the general solution of the 
dissipation inequality in the form of dynamic equations |22| . 

Here, in this paper we have chosen more direct approach with assuming suitable specific forms of the entropy l|II.4|l 
and the entropy current (|II.6p . The resulted entropy production makes possible to derive explicit evolution equations 
for the fluxes. The procedure allows to close the system and this is due to the Second Law. 

The expressions (IIIL2|I and (|IIL3p represent a system of 12 differential equations which, together with the balance 



of energy (jll.ip . allows - in principle - the determination of the 13 unknown functions e, qi and ^ij. However it is not 
in the balance form and hence, to achieve that task, some additional assumptions seem to be necessary. This subject 
will be investigated next. 



IV. BALANCE FORM FOR EVOLUTIONARY SYSTEMS 



Up to now we did not restrict the inductivity tensors m and n at all. Let us assume now, that they are constant 
and non-degenerate. Let us remark, that non-degeneracy is not trivial, since there exist real materials for which it 
is not guaranteed. A classical example is given in |23j . where an electric circuit described by dynamic variables is 
considered. 

For the sake of simplicity, let us assume that the inductivities take the form 

m = —ml, n — — nl, (IV. 1) 

where m and n are positive real coefficients and I means the unitary tensor in the corresponding tensorial space. 
Finalh^ according to arguments from the kinetic theory of rigid heat conductors, we assume that $ is symmetric 
P, 12, • Under the simplification above the equations lllll.2|) and 1111.3^ yield 

1 / d-^sn \ 

-q^~'^'^k'k = —[Likqk + Miki'^ki--^e,i], IV.2 
TO V oe-^ I 

1 / \ vn 

= -[Ptjkqk + Qijki'^kij + —qz'j- (IV.3) 

Here we denoted the symmetries of the constitutive functions, too. From (|IV.2p - ljIV.3|) results in that the evolution 
of q is ruled by a set of balance laws with gradient dependent source terms while tensor $ obeys a set of ordinary 
differential equations (with a given q). The system (|IL1|I . (jIV.2|) and (|IV.3|I may be regarded as the four moments 
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theory for rigid heat conductors of Extended Thermodynamics provided the equations ljIV.3|) are interpreted as 
suitable closure conditions 0,0. This point of view is quite different from that often applied in the current literature, 
where the closure relations are deduced either by different assumptions inside the kinetic theory by a maximum 

entropy formalism '25^ or by some algebraic considerations '2^ . To our opinion it is remarkable that closure relations 
have been obtained in a fully thermodynamic framework. 

Let us prove now that the 10 moments theory for rigid heat conductors may be obtained by a model with 27 
dynamic variables. Therefore in the following we assume that the introduced dynamic variables are symmetric. In 
this way we can compare our results with those of the kinetic theory. We assume the constitutive equation pi.3|l 
takes the form 



(IV.4) 



where ^[ij] is a symmetric second order tensor and "^[ijk] denotes a component of a fully symmetric third order tensor 
^ which enters the basic field. Here and in the following the bracket [jjj denotes a fully symmetric tensor in the 
bracketed indices. As a consequence we can write the entropy density as follows 



(IV.5) 



s{e,qi,^[ij],'^[ijk]) = so{e) + ^m^jq^qj + ^nijki^[tj]^[ki] + ^''jjfeipg*fofe]*[ip?], 

where the meaning of the negative definite tensor r is plausible. Here, in addition to the symmetry properties pL 
we can see, that riijki — njiki and riijki = riijik that we could denote as ?^[[ij][fc;]] ■ Similarly, the symmetry properties 
of r could be denoted by fHijmipq]] ■ In the following we omit the notation of the symmetries m, n and r. Moreover, 
if m, n and r are constant, we get 



dso 



ds 



-qt 



ds 



Due to (|1V.5|I and (|1V.6|I the entropy production takes the form 



*[feH] = -Q^qi + ruhk^if^i^qk + nijhk^ihk]^[iji], 



(IV.6) 



rrihpqp'k + nhkim'^ [imp]' p + nhkim^[im] 



[hk] 



nijhm^[h7n]'k + rijklpq^[lpq\ 



[ijk] ■ 



(IV.7) 



Again, due to (jlV.Sp the matrices m, n and r are symmetric (and hence invertible) and we are allowed to suppose 

r — — rl 



that m and n take the form (lIV.lll while 



with r positive and I unitary. Finally we get 



9% 



(IV.. 



(IV.9) 



-q{vj\ 



(IV.IO) 



(IV. 11) 



Here E, R, S, T and Z are Onsagerian conductivity tensors defined on the thermodynamic state space. In this way 
we get the 10 moments theory by (jlLlll . (I1V.9P and (jlV.10|l together with the closure relation (jlV.ll|l . It is worth 
noticing that the previous procedure of the extension of the basic fields can be continued analogously. That represents 
a general method to obtain extended thermodynamic systems together with suitable closure relations. 
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V. EXTENDED THERMODYNAMICS OF SECOND SOUND 

Second sound, i.e. thermal wave propagation, is a typical low temperature phenomenon which can be observed, for 
instance, in dielectric crystals such as Sodium Fluoride (NaF) and Bismuth (Bi) p^ . . Its appropriate description 
requires an extension of the classical Fourier's theory leading to the paradox of an infinite speed of propagation of 
thermal disturbances [T^ . From the microscopic point of view heat transport at low temperature is modelled through 
the phonon gas hydrodynamics [29l Isof. Phonons are quasi-particles which obey the Bose-Einstein statistics. In a 
solid crystal at low temperature they form a rarefied gas whose kinetic equation can be derived similarly to that of 
an ordinary gas. Phonons may interact among themselves and with lattice imperfections through two different types 
of processes: 

i) Normal-(N) processes, that conserve the phonon momentum; 

ii) Resistive-(R) processes, in which the phonon momentum is not conserved. 

The frequencies and vr of normal and resistive processes determine the characteristic relaxation times tm = 

and tr = Diffusive processes take over when there are many more R-processes than N-processes. If instead there 
are only few R-processes and many more N-processes, then a wave like energy transport may occur. This is said second 
sound propagation. Finally, there is a third mechanism of energy transport due to the presence of ballistic phonons 
travelling through the crystal without any interaction. The three different mechanisms of energy transport described 
above can be represented at a different level of approximation in the framework of Extended Thermodynamics. To 
achieve this goal let us suppose the dynamic variables are the first moment pi, {i = 1, 2, 3) and the momentum flux 
N[ij], = 1,2,3). By phonon gas hydrodynamics follows that the functions pi are connected to the heat flux 
by the relation qi = c^pi where c means the Debye's phonon velocity. Moreover, from kinetic theory we apply the 
interrelation of the traces of consecutive currents, which in our case is expressed by the definition N^^-^ := e. Therefore 
it is convenient to decompose the momentum flux N^^j^ into an isotropic part and a deviatoric part according to the 
equation 

A^M = ^e,5„ +7V<y>, (V.l) 

where Af<y > is symmetric and traceless. Due to the decomposition above the energy density coincides with the trace 
of N[ij] . As a consequence, we have only nine independent thermodynamic variables (9-field model) whose evolution 
is controlled by the balance laws 

P^+N[,kyk = P^, (V.2) 

N[ij] + J[ij]k'k = P[ij]- (V.3) 

Here Pi and P[y] are suitable production terms while J[ij]k represents the flux of iV[y], to be assigned through a 
closure relation. Under suitable constitutive assumptions on J[ij]k and P[ij], the balance of energy 

e + c'p^r, = (V.4) 

can be obtained by taking the trace of ljV.3p . Later on we illustrate such a property in more details by reinspecting 
two typical examples often used in the applications. 

4-fleld model 

This model holds true under the approximation tn = and describes second sound effects without taking care for 
N-processes. It is controlled by the equations 

e + c^pe^ = 0, (V.5) 



Pi + 7:e'i = Pi, (V.6) 

3 Tfi 



If one assumes <I>[jj] — c^eSij then the system above is recovered by ljIV.2p and l|IV.3|l . Moreover, if in ljIV.2p we 
assume the isotropy of the constitutive equations 

Mm=0, L,j = —6,j, (V.7) 

TR 
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e'l = Pi, (V.8) 



then we get 

On the other hand, due to the high speed of phonons (^42 x W^cmsec'^ in NaF 'sflj, the constant mc^ is very 

big. Hence, the coefficient in IjV.Sp can be neglected whenever the absolute temperature T — (^ff^ has 

no jumps, i.e. whenever the Lax conditions for shock wave formation [sH Is^ ] are not fulfilled. In such a case l|V.8|) 
reduces exactly to (jV.6|) . Further assuming isotropy and some more we get 

P[xj]k = 0, Q[ij][ki] =0, — = (V.9) 

and taking into account the decomposition P[iij] = ^Pk'k^ij + P<i'j>, with p<i'j> symmetric and traceless, equation 
(|IV.3p can be rewritten as follows 

^eSij + ^Pk'kSij + c^p<i'j> = 0. (V.IO) 
Finally, taking the trace of (|V.10|I we get the balance of energy (|V.5|I . 
9-field model 

Such an approximation holds for small tn and arbitrary r^. It is able to describe second sound taking into account 
the effects due to the propagation of ballistic phonons. Again, the calculated values of wave speeds are qualitatively 
but not quantitatively in agreement with the measured ones. The resulting system of equations is 

e + cV. = 0, (V.U) 



K + ^e-, + iV<,y>,, = p„ (V.12) 

O TR 

N<^3> + c^\p<vj> = --iV<y>, (V.13) 
o r 

where 

+ ^ (V.14) 

T TR TN 

is the total collision frequency. Let us observe that equation (|V.13p is not in the balance form since there is no flux 
of iV<ij> inside. Hence, nothing prevents to regard it as a closure relation and use ljIV.2|) and (|IV.3p to derive the 
governing system of equations. Then, let us make the identification 

= ye(5ij + c27V<y>, (V.15) 

and let us assume isotropic constitutive equations 

TTl 1 1 777. 

Mi[ki] = 0, Lij = — Sij, P[,j]k = 0, Q[ij][ki] ^ — ^ijSki + —SikSji — = c^. (V.16) 

TR T T Tl 

Here Q[ij][ki] is a constant isotropic tensor. There are only two material quantities because of the symmetries of 
Q[ij][ki]- Here we require that ^ = ^. Then by pV.2|l and ljIV.3(l we get 



Px 



V3 mc' ae^ / tr 
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1 • i 2 

-eOij + A^<y> + -pk'kOij + c p<,i'j> = N<:ijy. (V.18) 

O O T 

The trace of (jV.lSp yields the balance of the energy (jV.lip . As a consequence, due to (|V.11|) the last equation reduces 
to 

N<^J> + C^P<^'J> = --A^<y>. (V.19) 

r 

Also in this case in the absence of shocks (|V.17ll reduces to (|V.12|I while l)V.19|l differs from l)V.13|l only for the 
coefficient in front of p<i/j>. Such a discrepancy is of no concern at that level since the 9-field model gives only 
qualitative agreement with the experiments 



VI. CONCLUSIONS 



We reinspected the problem of thermal wave propagation and proved that the balance structure can be recovered 
with the Verhas form of the entropy current (|II.7|) if the thermodynamic inductivity coefficients and riiji^i are 
constant in the general entropy function (|II.4|I . The results above prove that the balance form of evolution equations 
is compatible with the mathematical structure of classical IT. The balance structure was obtained only for a subset 
of the independent thermodynamic variables. In fact, the highest order fluxes are controlled by ordinary differential 
equations, here regarded as closure relations. 

For our further investigations we have assumed some interrelations of the consecutive currents IjV.lj) , to be com- 
patible to the moment interpretation of the kinetic theory. With assuming isotropic material we have arrived to the 
expected system of equations of the 4 and 9 moment theory of rigid heat conductors, thereby demonstrating the 
compatibility of the two approaches. Let us observe the minor, but important additional assumptions (e.g. r = r'). 

Nowadays in ET the balance form of the evolution equations is a well motivated assumption from the kinetic 
theory and the entropy current is a derived quantity with a restrictive potential structure (|l.l|l . Moreover, the derived 
restrictions can be exploited predictively in different particular cases (e.g. for monatomic gases p. 52-60. We do 
not know related calculations for rigid heat conductors) . In ET one of the important problems is to give a reasonable 
closure of the arising hierarchy. 

However, the applications of the phenomenological theory in solids and fluids (see e.g. Q) and recent generalizations 
of the principles to incorporate nonlocality are far beyond the validity of the kinetic theory. That was our most 
important motivation of the recent investigations. In this paper we have shown that with suitable constitutive 
assumptions on the form of the entropy current and on the entropy one can recover the balance structure. We have 
demonstrated that one can give a kind of natural closure of the hierarchy. 

On the other hand, the conditions of using all the currents as state variables and to restrict ourselves to constant 
inductivities looks like too restrictive. However, to assume balance form evolution equations, where the current is a 
constitutive function is a weakly nonlocal theory from our point of view. E.g. 

ii = -^ik.kie, qi) = /(e, g^, e>i, qej). 

Calculation of the derivatives in the divergence of $ gives that the form of the evolution equation, defined by the 
function /, depends on the space derivatives of the state variables, too. Therefore, for a systematic investigation of the 
interrelation of the evolution equations and the entropy current without imposing the balance structure one should 
enlarge the basic thermodynamic space. Hence, one should weaken an other basic assumption in ET, the locality, 
by including the gradients of the state variables in the constitutive state space [l^, ■ ^ study of the effects of 
nonlocality on the evolution of higher order fluxes is developed in [s^l • 

Simplifying the arguments that we wanted to collide here we can say that although from the point of view of kinetic 
theory the specific form of the entropy current is a restriction, from a phenomenological point of view the balance 
structure looks like too restrictive. 
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